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ABSTRACT 

<N 

» I , Aims. The aim is to investigate the star-formation and LINER (Low Ionization Nuclear Emission Line Region) activity 

within the central kiloparsec of the galaxy NGC 1614. In this paper the radio continuum morphology, which provides a 
tracer of both nuclear and star-formation activity, and the distribution and dynamics of the cold molecular and atomic 
gas feeding this activity, are studied. In particular, the nature of an R~300pc nuclear ring of star-formation and its 
ly-^ , relationship to the LINER activity in NGC 1614 is addressed. 

Methods. A high angular resolution, multi-wavelength study of the LINER galaxy NGC 1614 has been performed. Deep 
observations of the CO 1-0 spectral line were performed using the Owens Valley Radio Observatory (OVRO). These 
data have been complemented by extensive multi-frequency radio continuum and Hi absorption observations using the 
Very Large Array (VLA) and Multi-Element Radio Linked Interferometer Network (MERLIN). 

Results. Toward the center of NGC 1614, we have detected a ring of radio continuum emission with a radius of 300 pc. 
This ring is coincident with previous radio and Paa observations. The dynamical mass of the ring based on Hi absorption 
is 3.1 x 10 9 Mq. The peak of the integrated CO 1-0 emission is shifted by 1" to the north-west of the ring center. An 
upper limit to the molecular gas mass in the ring region is ~ 1.7 x 10 9 Mq. Inside the ring, there is a north to south 
elongated 1.4 GHz radio continuum feature, with a nuclear peak. This peak is also seen in the 5 GHz radio continuum 
and in the CO. 

Conclusions. We suggest that the R = 300 pc star forming ring represents the radius of a dynamical resonance - as 
an alternative to the scenario that the starburst is propagating outwards from the center into a molecular ring. The 
ring-like appearance is probably part of a spiral structure. Substantial amounts of molecular gas have passed the radius 
of the ring and reached the nuclear region. The nuclear peak seen in 5 GHz radio continuum and CO is likely related 
to previous star formation, where all molecular gas was not consumed. The LINER-like optical spectrum observed in 
NGC 1614 may be due to nuclear starburst activity, and not to an Active Galactic Nucleus (AGN). Although the 
. presence of an AGN cannot be excluded. 

Key words, galaxies: evolution — galaxies: individual: NGC 1614 — galaxies: starburst — galaxies: active — radio 
Q\ ■ lines: ISM — ISM: molecules — ISM: atoms 
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O ' !■ Introduction Previous studies of NGC 1614 include e.g. an optical, 

<^> near-infrared, radio continuum and hi study by Neff et 

. NGC 1614 is a luminous infrared galaxy (LIRG) at a dis- a l. (1990). They suggested that the spectacular structure 

| tance of 64 Mpc (for Hq=75 kms _1 Mpc _1 ). The galaxy is with tidal tails or plumes is the result of an earlier inter- 

• *"H , barred and interacting (morphologically classified as type action with another galaxy of comparatively modest mass 

' SB(s)c pec) and its nuclear optical spectrum shows both and impact velocity. Alonso-Herrero et al. (2001), hereafter 

H | starburst and LINER (Low Ionization Nuclear Emission AH2001, have investigated NGC 1614, in the optical and 

Line Region) activity. Although there is still no single con- near-infrared. They detected a starburst nucleus of about 

sensus to what is powering the emission in LINERs, nu- 45 pc in diameter based on deep CO stellar absorption, 

clear starburst or AGN (Active Galactic Nucleus) activity The nucleus is surrounded by a ^600 pc diameter ring 

are two major candidates (e.g. Ho 1999; Terashima et al. of current star formation, which is revealed in Paa line 

2000; Alonso-Herrero et al. 2000). Bar-driven inflow of gas emission. Just outside the star forming ring, a dust ring 

in galaxies has been suggested to trigger nuclear starbursts is indicated by its extinction shadow in H — K. Neither 

as well as feeding AGNs (e.g. Simkin, Su & Schwarz 1980; Neff et al. (1990) nor AH2001 find any indication of an 

Scoville et al. 1985). However, the radial inflow of gas along AGN. Scoville et al. (1989) used the Owens Valley Radio 

the bar may be slowed down at certain radii, often associ- Observatory (OVRO), which contained three elements at 

ated with inner Lindblad resonances (ILR) (Combes 1988; that time, to map NGC 1614 in CO 1-0 at a resolution of 

Shlosman, Frank & Begelman 1989). 4" x 6". They found an unresolved CO-concentration with 

a mass of 6 x 10 9 M . 

We have used OVRO to study the distribution and dy- 

Send offprint requests to: S. Aalto namics of the molecular gas at 2" resolution. We have also 
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Fig. 1. Overlay of OVRO CO 1-0 integrated intensity contours over a F606W WFPC2 HST image. The contours are in 
percent of the peak value of 34 Jykms" 1 per beam of 2."75 x 2."40. The levels are 8, 11, 16, 23, 32, 45, 64 and 90%. 
The greyscale of the optical image is logarithmic and arbitrary. The astrometrical alignment of the HST image with the 
CO map is within l."5 (see Sect. 13. 1|) . 



studied the morphology of the 1.4 GHz radio continuum 
as well as the distribution and dynamics of the neutral gas 
via Hi in absorption at arcsecond resolution, using the Very 
Large Array (VLA). To study the morphology of the nu- 
clear radio continuum, we used the Multi Element Radio 
Linked Interferometer Network (MERLIN) and obtained 
maps at higher (0.5 arcsecond) resolution, at 1.4 GHz and 
5 GHz. The purpose of the studies was to provide informa- 
tion about the feeding and nature of the central activity 
giving rise to the LINER like spectrum in NGC 1614. In 
particular, we would like to address the question whether 
the LINER activity is due to an AGN or to a nuclear star- 
burst. 

2. Observations and data reduction 

2.1. Owens Valley Radio Observatory 

The OVRO mm interferometer was used to map NGC 1614 
in CO 1-0 in the inner 60" (R — 9.3 kpc). The array con- 
sists of six 10.4 m telescopes and was used in the equatorial 
and high resolution configuration. These observations were 
carried out in February 1996. The NRAO's MPS software 
package was used to deconvolve the images. We used data 
from both configurations together to produce maps with 
high sensitivity and image fidelity, and data from the high 



resolution configuration only to produce maps with high 
angular resolution. The resulting low and high resolution 
synthesized beams were 4. "44 x 4. "10 and 2. "55 x 2. "40. 
The primary beam diameter was 60". The digital corre- 
lator was centered at 113.4 GHz (4800 kms~ ) and was 
configured to cover 448 MHz (1200 krns -1 ), with 8 MHz 
(22 krns -1 ) resolution. Typical system temperatures were 
400 K. 

2.2. Very Large Array 

NGC 1614 was observed at 1.4 GHz with the VLA in the 
A-configuration in April 2006. The observations were inter- 
spersed with regular observations of the nearby phase cali- 
brator 0423-013. The flux calibrator was 3C 286. The obser- 
vations were centered on 1398 MHz, which corresponds to 
a velocity of 4778 krns -1 . A two-IF mode was used, with 
a band width of 6.25 MHz per IF. The IFs were parallel 
in frequency, which resulted in a total frequency coverage 
of 6.25 MHz (1360 krns" 1 ). Each IF consisted of 32 chan- 
nels with a bandwidth of 96 kHz, which corresponds to a 
velocity resolution of 21.3 krns -1 . The data were edited 
and calibrated with AIPS standard procedures. The cali- 
brated data were Fourier transformed to produce spectral 
line data cubes with two different weightings. The line free 
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Fig. 2. OVRO CO 1-0 integrated intensity, velocity field and major axis (north to south) position velocity diagram. The 
integrated intensity (left panel) contours are in percent of the peak value of 34 Jykms -1 per beam of 2. "75 x 2. "40. 
The levels are 8, 11, 16, 23, 32, 45, 64 and 90%. The greyscale range is from to 34 Jykms -1 beam -1 . The velocity 
field (middle panel) greyscale and contour range is from 4700 to 4850 kms -1 , with 25 kms -1 contour increments. The 
position velocity diagram (right panel) is averaged over a 2. "5 slit, centered on the peak of the integrated intensity (see 
Table 1). The dotted line represents the rotation curve in this region, and the arrow points at an additional kinematical 
component. The position velocity diagram contours are 40, 50, 60, 70, 80 and 90% of the peak value of 0.28 Jy per beam. 
The dotted ring marks the position of the starburst ring. Note that it is offset to the east from the CO peak by 1". 



channels were combined to form continuum images, which 
were then subtracted from the spectral line data. The con- 
tinuum images and the continuum subtracted spectral line 
data were separately deconvolved with the AIPS task AP- 
CLN, and were then recombined to form the deconvolved 
spectral line data cubes. Two different weighting schemes 
were applied to the data in order to obtain maximum sensi- 
tivity and angular resolution. For maximum sensitivity, nat- 
ural weighting was applied, which resulted in a synthesized 
beam of 2. "33 x l."42. For maximum angular resolution, 
a more uniform weighting (robustness parameter —2) was 
applied. The resulting synthesized beam was l."39 x l."01. 
Further analysis such as presenting cleaned contoured con- 
tinuum images, absorption spectra and moment maps were 
done with standard AIPS tasks. 

In order to provide complementary high resolution ra- 
dio continuum imaging of the central region of NGC 1614, 
data at 5 and 8.4 GHz was obtained from the VLA archive 
and re-imaged. These short snapshot data were observed 
on 27th July 1999 in the VLA's highest resolution A- 
configuration. Each data-set was was calibrated using stan- 
dard data reduction techniques, including phase referencing 
using a nearby calibrator, within the AIPS packages. The 
flux density scale of these observations was calibrated with 
respect to 3C286 using the Baars et al. (1977) scale. 

2.3. Multi Element Radio Linked Interferometer 

MERLIN was used in August 2005 to observe NGC 1614 
at 1420 MHz. The observations were interspersed with reg- 
ular observations of the nearby phase calibrator 0436-089. 
3C 286 was used as the primary calibrator and 0552+398 
as the secondary. They were both observed at the begin- 
ning and end of the observing run. Dual bands of circu- 
lar polarization were recorded over a total bandwidth of 8 
MHz, which was correlated into 64 channels width a band- 
width of 125 kHz each, which equals a velocity resolution 
of 26 kms -1 . Initial editing and calibration of the data 
was done at Jodrell Bank in September 2005 using the lo- 



cal MERLIN DPROGS software. These data were read into 
AIPS and further calibration was done using the MERLIN 
pipeline, which included several cycles of self calibration 
on the phase calibrator. Our target source was not itself 
suitable for self calibration, following the constraints given 
in the MERLIN handbook, so the phase corrections de- 
rived from the phase calibrator were applied to our target 
source in the MERLIN pipeline. The calibrated uv data- 
set was Fourier transformed with no deconvolution initially 
applied. The line free channels were combined to produce 
continuum images which were used to subtract the con- 
tinuum contribution in the spectral line cubes. The con- 
tinuum images and the continuum subtracted spectral line 
data cubes were separately deconvolved with the AIPS task 
APCLN, and were then recombined to form the deconvolved 
line data. The rather low declination of this source resulted 
in an elongated synthesized beam of 0."50 x 0."17. 

In addition to the 1420 MHz MERLIN data, MERLIN 
observations at 4994 MHz were also obtained. These data 
consist of two full track observations of NGC 1614, made on 
24th April and 19th May 2000. In each case observations of 
NGC 1614 were interspersed with scans of the phase refer- 
ence source 0436-089, with additional observations of stan- 
dard point source and flux density calibration sources at 
either end of the observing run. In total NGC 1614 was ob- 
served for 12 hr over the two observing dates. Both of these 
observations were independently calibrated and imaged us- 
ing standard routines before the data were combined for 
final imaging of the target. 

MERLIN has a shortest baseline spacing of ~11 km. As 
a consequence 5 GHz observations with MERLIN are in- 
sensitive to diffuse radio structures larger than ~1 arcsec. 
In order to restore these missing short spacing data, our 
MERLIN 5 GHz data were combined within the uv plane 
with the calibrated VLA A-configuration snapshot data 
(see Sect. 12. 2p . applying appropriate weightings to each 
data set according to their sensitivity. The resulting com- 
bined data set was then deconvolved and imaged to provide 
maps both sensitive to the large diffuse radio emission and 
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Table 1. Adopted properties of NGC 1614. 
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Systemic velocity 


4778 kms -1 


Distance 13 


64 Mpc 


Spatial scale 


1" = 310 pc 


Inclination c 


51° 


r d 

Mr 


3 x 10 11 L Q 


Adopted conversion factor 


2.3 x 10 20 cm" 2 (K kms -1 ) -1 


Adopted velocity convention 


Optical heliocentric 



a) : Neff et al., 1990 

b) : For #o=75 kms -1 Mpc -1 

c) : For the R=l" Pa ring, Alonso-Herrero et al., 2001 

d) : Alonso-Herrero et al., 2001 



with an angular resolution intermediate between the indi- 
vidual MERLIN and VLA images. 
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3. Results 

3.1. The OVRO CO 

Figure [T] shows the high resolution integrated intensity 
contours overlayed on an F606W WFPC2 (optical) HST 
(Hubble Space Telescope) image. The optical peak posi- 
tion has a positional uncertainty of less than 1" (Neff et 
al., 1990) and we estimate that the astrometrical alignment 
between the optical and CO is better than l."5. The same 
CO integrated intensity map is also shown in Fig. [2] (left 
panel). The CO is centrally peaked, with an extension to 
the northeast, with CO at lower integrated intensity. A sec- 
ondary peak is detected at ~ 3" (1 kpc) to the northeast. 

The maximum integrated intensity is 34 Jy kms - , in a 
beam of 2. "75 x 2. "40. With a standard CO to H 2 conversion 
factor (see table [1]) , this corresponds to a molecular mass 
of 1.3 x 10 9 M Q in the central beam. The projected surface 
density is 2600 M Q pc -2 in the central beam. The depro- 
jected surface density is ~ 1600 M pc -2 for the adopted 
inclination of 51°. 

The high resolution velocity field map is shown in Fig. [2] 
(middle panel). The velocity contours are perpendicular to 
the main axis of the central CO structure, with velocities 
ranging from 4675 to 4825 kms -1 . On a linear scale, the 
observed velocity gradient in the central 2" corresponds 
to ~ 140 kms -1 per kpc. In the extension to the northeast, 
the velocity contours are more irregular, probably due to an 
additional kinematical component in this region. In Fig. [TJ 
this feature is spatially consistent with a crossing dust lane. 

The position velocity diagram, along the north to south 
direction, is shown in Fig. [5] (right panel). The position 
velocity diagram is averaged over a 2. "5 slit, centered on the 
peak of the integrated intensity. Three peaks are resolved. 
The two strongest peaks occur to the north and south of 
the center. The third peak, which is weaker, occur toward 
the center of the CO integrated intensity. 

The low resolution CO data do not show any additional 
features as compared to the high resolution maps, but have 



Fig. 3. Robustly weighted VLA radio continuum at 1.4 
GHz. The synthesized beam is l."39 x l."01. The noise 
level is 0.32 mJy beam -1 . The contours start at 1 mJy and 
increase with a factor of V2 per level. The greyscale range 
is from to 25 mJy beam -1 . 



higher sensitivity. The total flux recovered is 80 Jy kms - , 
which corresponds to a molecular mass of 3.0 x 10 9 M Q . 
NGC 1614 was previously observed with OVRO (Scoville 
et al., 1989) at lower resolution (4" x 6"), where a total of 
103 Jy kms~ was detected. 

The global star formation rate (SFR) is related to Lfir 
(£fir = 3 x 10 11 L Q ) and is given by Kennicutt (1998) as 
SFR=Lfir/5.8 x 10 9 M Q yr -1 . This relation indicates an 
SFR of - 52 M Q yr -1 in NGC 1614. The star formation 
efficiency, SFE, is given by SFE=£ F ir/M H2 an ^ is ~ 100 
(in the unit Lq/Mq ). This is in good agreement with the 
expected value of ~ 100 for interacting galaxies, given by 
e.g. Young et al. (1989); Planesas et al. (1997). 

3.2. The VLA and MERLIN Radio Continuum 

Figure [3] shows the VLA robustly weighted 1.4 GHz radio 
continuum map. The synthesized beam is l."39 x l."01. 
The total flux density recovered is 130 mJy and the noise 
level is 0.32 mJy beam -1 . The maximum flux in one beam 
is 29 mJy. In a naturally weighted 1.4 GHz radio continuum 
map (with a synthesized beam of 2. "33 x l."42), we recover 
a total flux density of 140 mJy. 

Figure H shows the MERLIN 1.4 GHz and 5 GHz radio 
continuum of the central region of NGC 1614, as well as 
archival VLA 5 GHz and 8.4 GHz radio continuum of the 
same region. Figure [S] shows the MERLIN 5 GHz radio 
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Fig. 4. Upper left panel: MERLIN 1.4 GHz with a synthesized circular beam of 0."3. The greyscale is from to 4 mJy 
beam -1 . The first contours are at —0.32 and 0.32 mJy beam -1 and then increase with a factor of v2 per level. Upper 
right panel: Combined MERLIN and VLA 5 GHz with a synthesized beam of 0."34 x 0."28. The greyscale is from to 
2.5 mJy beam -1 . The first contours are at —0.1 and 0.1 mJy beam -1 and then increase with a factor of y/2 per level. 
Lower left panel: VLA 8.4 GHz with a synthesized beam of 0."34 x 0."23. The greyscale is from to 1.4 mJy beam -1 . 
The first contours are at —0.029 and 0.029 mJy beam -1 and then increase with a factor of \/2 per level. Lower right 
panel, greyscale: MERLIN 1.4 GHz with a synthesized circular beam of 0."3. The greyscale is from to 4 mJy beam -1 . 
Contours: MERLIN 5 GHz with a synthesized circular beam of 0."2. The first contours are at —0.32 and 0.32 mJy 
beam -1 and then increase with a factor of y2 per level. 



continuum overlayed on a Paa mapQ made with archival 
data from the Hubble Space Telescope, HST.. 

1 Based on observations made with the NASA/ESA Hubble 
Space Telescope, obtained from the data archive at the 
Space Telescope Science Institute. STScI is operated by the 
Association of Universities for Research in Astronomy, Inc. un- 
der NASA contract NAS 5-26555 



The Paa (the greyscale in Fig. [5]) traces out a ring-like, 
patchy structure with a radius of ~ 1" (310 pc). This has 
previously been reported by AH2001. Inside of this ring, the 
Paa intensity is lower, except at the center, where there is 
a Paa peak. We have spatially aligned this peak with the 
position of the central peak of the 5 GHz MERLIN radio 
continuum (the contours in Fig. [5]). 
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Fig. 5. Overlay of MERLIN 5 GHz in contours over a 
HST Pact map in greyscale. The synthesized beam of the 
MERLIN data is circular with a diameter of 0."2. The 
dashed contour is at —0.32 mJy beam -1 , and the positive 
contours start at 0.32 mJy beam -1 and has a contour incre- 
ment of a factor of \/2 per level. The greyscale is arbitrary. 



The MERLIN 5 GHz map, the combined MERLIN and 
VLA 5 GHz map and the VLA 8.4 GHz map all trace out 
the same ring-like structure as the Paa (see Fig. 2]). We have 
chosen to refer to this structure as the star forming ring. 
The combined MERLIN and VLA 5 GHz map (upper right 
panel, Fig. 2]) coincides remarkably well in space with the 
Paa, around the whole circumference of the star forming 
ring. In the MERLIN map alone (Fig. [5] and lower right 
panel in Fig. 2]), some of the emission is resolved out, and 
mainly the eastern and western parts of the star forming 
ring are detected. It is also detected in the VLA 8.4 GHz 
map (lower left panel, Fig. 2]). 

The MERLIN 1.4 GHz map (Fig-H upper left panel) re- 
veals a north to south bar like structure, with several peaks 
resolved. The two brightest peaks occur 1" (310 pc) from 
the center, on the star forming ring toward the north and 
south respectively. The southern peak is spatially extended 
to the north, toward the center. The third peak occurs to- 
ward the west, in the Paa ring. 

There is an unresolved peak at the center of the star 
forming ring, visible in most of our radio continuum maps. 
In particular, the peak is isolated in the MERLIN 5 GHz 
map, and tentatively present in all maps. 

3.3. The VLA Hi Absorption 

We have detected Hi in absorption in the central R w 2" 
of NGC 1614. Figure [6] (lower panel) shows spectra from 
the VLA continuum subtracted spectral line data cube. Hi 
spectra from the MERLIN data (not shown) are consistent 
in width and optical depth with the VLA spectra, but are 
not used for the analysis or moment maps due to low signal 
to noise ratio. 



The separation between individual spectra in Fig. [5] is 
0."5 in both RA and dec. The absolute depth of the ab- 
sorption does not only depend on the column density of 
foreground Hi, but also on the flux density of the back- 
ground radio continuum. Hence the absorption is deeper in 
the central spectra, and no absorption is detected outside 
a radius of l."5 in the east to west direction. To the north 
and south (not shown in this map) there are tentative de- 
tections out to a radius of 2". 

The moment maps of the Hi data are shown in Fig. [7J 
The left panel is the integrated absorption. Its peak oc- 
curs close to the radio continuum peak ( but is not a true 
measure of the Hi distribution). The middle panel shows 
the integrated opacity (J rdV). This map has been blanked 
outside of the region where the radio continuum flux den- 
sity is less than 10 mJy beam -1 , since the noise in this 
map increases dramatically where there is little or no ra- 
dio continuum. The integrated opacity is a measure of the 
foreground absorbing Hi, and there is a clear east to west 
gradient in this map. 

The column density of absorbing atomic hydrogen, Nn, 
has been calculated at several positions using the equation 

7V H = 1.823 x 10 18 T spin / rdV atoms cm -2 

where T sp i n is the spin (excitation) temperature, and J rdV 
is the integrated opacity (from Fig. [7]) in kms -1 . T sp i n can 
be reasonably assumed to be 100 K (Maloney, Hollenbach 
& Thielens, 1996), although it should be noted that the 
value of T S pi n is dependent on the physical conditions of 
the gas. 

The maximum column density occur to the west and 
northwest of the radio continuum peak, where the peak 
opacity is ~ 0.25 and the velocity width is ~ 200 kms -1 . 
These values result in a column density approaching 
10 22 atoms cm -2 . East of the radio continuum peak, the 
column density is approximately 50% lower. 

The right panel of Fig. [7] shows the velocity field. There 
is a clear north to south velocity gradient, where the south- 
ern end is approaching and the northern end receeding, 
similar to what is seen in CO. 

3.4. Dynamical results with GAL 

We used the velocity of the OVRO CO data to fit a solid 
body rotation curve to the inner R — 5" with the MPS task 
GAL. For solid body rotation, the inclination and rotational 
velocity can not both be fitted with GAL, since they are 
directly dependent on each other. We therefore kept the 
inclination fixed at i = 51°, which is the inclination of the 
star forming ring (AH2001). The dynamical center was kept 
fixed at the CO-peak (see Table [1]). The fit gave a rotation 
of 40 kms - per arcsecond (130 kms - per kpc), which 
we used to calculate a dynamical mass within a radius of 
4" (1240 pc). We used the Keplerian relation 

M dyn = 2.3x 1O 8 (^) 2 ( I | J )M 

where V ro t is in kms -1 and R in pc. This resulted in a dy- 
namical mass of 7.3 x 10 9 Mq, which has to be taken as 
an estimate since we can not rule out the presence of non- 
circular orbits. The ratio of dynamical mass over molecular 
mass detected in the integrated intensity map is ~ 2.5. 

The residual velocity field, after the fitted rotation curve 
was subtracted, was also calculated with GAL. The resid- 
ual velocities range from —40 kms -1 to 50km s - , and the 
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Fig. 7. The left panel shows the integrated VLA Hi absorption, derived from the continuum subtracted spectral line 
cube. The greyscale and contour range are from —500 to —1700 mJy bm _1 kms _1 . The contour increment is 200 mJy 
bnr'kms -1 . The synthesized beam is 2. "33 x l."42 (in all three panels). The middle panel shows the integrated optical 
depth (J rdV). The greyscale range is from to 65000 (in the unit kms~ ). The contour range is from 40000 to 60000 
with increments of 5000. Note that the errors in this map increase in areas where the radio continuum is weak, i.e. toward 
the edges of the map. To minimize this effect, this map has been blanked where the radio continuum is less than 10 mJy 
bm . The right panel shows the velocity field, derived from the continuum subtracted spectral line cube. The greyscale 
range and contours are from 4675 to 4850 kms _1 with 25 kms _1 contour increments. 



zero point in the residual field occur toward the dynamical 
center. 

The velocity field derived from the naturally weighted 
VLA continuum subtracted spectral line data cube is shown 
in Fig. [7] (right panel). The main north to south rotation is 
clearly seen, but also a deviation in the velocity field in the 
northern half of the galaxy. As seen in Fig. [5] (lower panel) , 
there are broad absorption line wings in several spectra. In 
order to isolate the main rotation from additional kinemat- 
ical components, we used the MPS task GAL to fit a solid 
body rotation curve to only the deepest, narrow absorption 
components. We used a cut off level of —3 mJy per beam 
to derive a main velocity field. We fitted a rotation curve, 
which was then subtracted from the original velocity field 
(Fig. El right panel). The fitted solid body rotation was ~ 
210 kms -1 per kpc, and the residual of the velocity field is 
shown in Fig. [6] (upper panel). The observed residuals are 
all blueshifted, with ~ 100 kms~ over most of the fitted 
area (the central R = 2"). 

4. Discussion 

4.1. Large scale molecular distribution and dynamics 

The integrated intensity and velocity field of the OVRO 
CO data is shown in Fig. [5] (left and middle panel). To 
the northeast, there is a CO extension, and the gradient in 
the velocity field falls, as compared to the gradient in the 
center. We believe that this extension is caused by a sep- 
arate, foreground gas component, associated with a dust 
lane, which can be seen in Fig. [TJ (the lighter region north- 
east of the nucleus). This component is likely not to be in 
the same plane as the spiral arms visible in the same figure. 
For further discussion of this feature see Sect. 14.81 

The position velocity diagram of the same data (Fig. [5J 
right panel) also shows an additional, blueshifted compo- 
nent which is associated with the foreground dust lane. The 
dotted line in the position velocity diagram was interpo- 
lated between the southern and central peaks (not contam- 



inated by the dust lane) and then extrapolated to the north. 
The additional component is marked in the figure. We con- 
clude that the bulk of the CO is associated with a sym- 
metric, slightly elongated bar like structure, and that there 
is a second CO component associated with a foreground 
dust lane to the northeast. The CO-contours in Fig. Q] also 
shows an extension to the west of the main CO peak. This 
extension is consistent with also being associated with the 
same foreground dust lane, which crosses the galaxy from 
the northeast to the southwest. 

Note that the CO peak is offset from the radio contin- 
uum peak (and also the starburst ring) by 1" to the north- 
west (Fig. O. This is consistent both with the east- west 
extinction seen by AH2001 and the east-west gradient in 
opacity in atomic hydrogen in absorption found in this pa- 
per in Sect. 13.31 This suggests that the extinction shadow 
found by AH2001 is not a part of a complete ring but actu- 
ally represents a true lopsidedness in gas and dust around 
the center. 

4.2. The atomic gas - column densities and dynamics 

As described in Sect. 13.31 and seen in the integrated opacity 
map (Fig. [71 middle panel), the maximum column density 
of absorbing Hi does not occur toward the radio contin- 
uum peak. To the contrary, the highest opacity occur at 
the edges of the integrated opacity map. This effect is bet- 
ter seen in the major axis position velocity diagram of the 
opacity (Fig. [5]), which clearly shows two separate opacity 
peaks, to the north and south respectively. The position ve- 
locity diagram is blanked where the background VLA radio 
continuum flux density is less than 10 mJy beam -1 , which 
spatially corresponds to ~ 500 pc to the north and south of 
the radio continuum peak. The peak opacities occur close 
to the edges of the position velocity diagram, i.e. at R « 
400 pc. This is consistent with, or slightly larger than, the 
radius of the star forming ring described in Sect. l3.2l and we 
conclude that most of the absorbing Hi is associated with 
the outer parts of this ring. The northern and southern 
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Fig. 6. The lower panel shows spectra of the VLA naturally 
weighted continuum subtracted spectral line cube. The cen- 
tral spectra are at the radio continuum peak (see Table [1} , 
and the separation between spectra is 0."5 in RA and DEC. 
The velocity scale is 200 kms -1 between ticks, and the tick 
marked with an arrow in the central spectrum is at 4800 
kms -1 . The absorption scale is from —8 to 2 mJy beam -1 , 
with increments between ticks of 2 mJy beam -1 . The up- 
per panel is the residual velocity field, after subtraction of a 
solid body rotation fitted only to the main (deeper than —3 
mJy per beam) kinematical component. The greyscale and 
contour range is from —110 to 10 kms _1 and the contour 
increment is 10 kms -1 . 



peaks are also seen in the integrated opacity map (Fig. 
middle panel), which also shows a clear east to west opac- 
ity gradient, with a peak to the northwest. This peak is 
spatially consistent with the western end of the dust lane 
(the light region northwest of the CO-peak in Fig. [TJ and 
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Fig. 8. Major axis (north to south) position velocity dia- 
gram of the VLA Hi opacity, r. The position velocity di- 
agram is averaged over a l."75 slit, centered on the radio 
continuum peak (see Table 1). The greyscale range is from 
0.075 to 0.2. The contour range is from 0.075 to 0.175 with 
increments of 0.025. 



is consistent with the extension seen in the CO-contours in 
this region. 

In Sect. 13.41 we isolated the main Hi rotation, and 
showed that the residual Hi velocity field is all blueshiftcd 
(Fig. with velocities of ~ 100 km s _1 over most of the fit- 
ted area (the central R = 2"). We interpret this as a large 
outflow, or super wind, associated with the past and current 
star formation. Superwinds are found in almost all infrared- 
luminous galaxies (Rupke et al. 2005). For an inclination of 
i = 51°, our observations correspond to a deprojected wind 
velocity of ~ 160 kms -1 . This is in good agreement with 
the average outflow velocity of the ionized gas in LIRGs of 
170 kms _1 (Lehnert & Heckman, 1996). 

4.3. The central radio continuum - a nuclear ring or spiral 



In Sect. 13.21 we described a star forming ring with a radius 
of 1" (310 pc), seen in the MERLIN and VLA radio con- 
tinuum and in the Paa archival data. (Fig. @] & It has 
been detected previously in the Paa, by AH2001, and with 
the VLA at lower resolution by Hummel et al. (1987) and 
Neff et al. (1990). 

The molecular mass in the central OVRO CO beam 
(2."75 x 2."40) is 1.3 x 10 9 M Q , using the conversion fac- 
tor in Table [TJ This mass can be used as an upper limit to 
the molecular mass associated with the star forming ring. 
This is a conservative upper limit since most of the CO 
peak appears not associated with the ring itself, but re- 
sides 1" to the north-west of it. To make a mass budget 
for the ring, we have estimated the dynamical mass in the 
central R = 1". The highest possible spatial resolution was 
needed, and we used the VLA Hi absorption data (with 
a synthesized beam of 2. "33 x l."42). Three peaks are re- 
solved in the position velocity diagram (Fig. [9]), and we 
estimated a projected velocity gradient in the ring of ~170 
kms _1 arcsec _1 , based on the slope of the dotted line in 
Fig. [9l A more conservative estimate would be to use the 
central line width at 50% of the maximum absorption in 
the position velocity diagram (the arrow in Fig. [H]), which 
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Fig. 9. Major axis (north to south) position velocity dia- 
gram of the VLA Hi absorption. The position velocity di- 
agram is derived from the continuum subtracted spectral 
line data cube and is averaged over a l."75 slit, centered 
on the radio continuum peak (see Table 1). The contours 
are at 40, 50, 60, 70, 80 and 90% of the peak absorption 
depth of —6.4 mJy per beam. The dotted line represents 
the rotation curve in this region. The arrow represents the 
central velocity width at 50% of the peak absorption depth. 

is 240 kms -1 , and would give a projected velocity gradient 
in the ring of ^120 kms _1 arcsec _1 . For these two veloc- 
ity gradients and an inclination of 51°, using the Keplerian 
relation. 

Af dyn = 2.3x 1O 8 (%) 2 ( I §)M 

where V TO t is in kms _1 and R in pc, we derive a dynamical 
mass range of 1.7 x 10 9 M Q (conservative) to 3.5 x 10 9 M 
in the central R = 1". 

This is a factor of ~ 1.5 — 2.5 higher than the value cal- 
culated by AH2001. They reported that the molecular gas 
mass, as calculated with a standard CO to H2 conversion 
factor, would clearly exceed the dynamical mass of the star 
forming ring. The higher angular resolution of our OVRO 
data, as compared to previous observations by Scoville et 
al. (1989), has better isolated the CO in the star form- 
ing ring, and provides a molecular gas mass estimate of 
1.3 x 1O 9 M0. This, and our higher dynamical mass, allows 
a mass budget including all of the molecular gas detected 
in the central OVRO beam. 

Note that the CO peak is offset from the radio contin- 
uum peak by about 1" suggesting that most of the molec- 
ular mass in the CO peak is not associated with the star 
forming ring. This further enforces the conclusion that there 
is no direct conflict between the estimated molecular mass 
and the dynamical mass of the ring. 

4.4. Is there a bar in NGC 1614? 

The origin of the star forming ring can either be a nuclear 
starburst, that has consumed most of the gas in the inner 
R = 1", and is progressing outwards, or it can have dynam- 
ical origin, possibly associated with an inner Lindblad res- 
onance (ILR), where the gas would pile up (Combes 1988a; 
Shlosman et al. 1989). NGC 1614 has previously been de- 
scribed as a textbook example of a propagating starburst 



that started in the nucleus of a late-type, large spiral galaxy, 
has grown outward to a radius of ~ 300 pc, and is poten- 
tially still growing into a circumnuclear ring of molecular 
material just outside this radius (AH2001). 

We suggest that, instead of being the result of a star- 
burst spreading outward from the center, the star forming 
ring represents the radius of a dynamical resonance. The 
ring-like appearance may be part of a spiral structure, ex- 
tending from the nuclear region, inside of the ring, out to 
the spiral arms visible in the F606W WFPC2 HST image 
(Fig. H]). 

This scenario is shown in Fig. [TU1 where the greyscale is 
the Paa emission and the different dynamical features are 
marked with arcs. The outer black arcs represent the main 
spiral arms, which are clearly visible in the F606W WFPC2 
HST image and Paa images (Fig. IB"! fe ITU]) . The dashed arcs 
represent the continuation of the main spiral arms toward 
the center, probably associated with the leading edge of a 
bar. We believe that the gas flows along the leading edge of 
the bar, and piles up at the radius of the star forming ring. 

Inside the ring, there may be a dynamically decoupled 
bar, in the north to south direction, which connects to the 
nucleus (the white arcs in Fig. [TUJ. The gas mass fraction 
in the central R = 1" is well over 10% which is the required 
lower limit (Hunt et al., 2008) for an inner bar to decou- 
ple from the outer bar. The north-south feature in the 1.4 
GHz MERLIN map (Fig. [JJ may represent such an inner, 
decoupled bar. 

The ring formed at the resonance acts as a barrier pre- 
venting further mass inflow, but Shlosman solved this by 
forming inner bars that may connect with the outer bar, 
i.e. nested bars (e.g. Shlosman 2003). A problem in our in- 
terpretation may then be that the nested bar forms after 
the ring formation and thus the activity in the center should 
therefore be younger than the activity in the ring. However, 
there are at least two arguments that counter this. Gas may 
flow to the nucleus before the ring is fully formed where it 
can feed an AGN or a starburst. The nuclear activity may 
not be triggered immediately but is dependent on gas build- 
up times and dynamics. Furthermore, the star formation in 
the ring itself may well be episodic and the ring could go 
through several bursts. 

Sarzi et al. (2007) discuss the formation of rings in sys- 
tems both with and without strong kpc-scale bars. The 
morphology of the resonant nuclear rings in their sample 
is almost identical to the ring of NGC 1614 complete with 
a nuclear AGN and/or starburst. Sarzi et al. discuss how 
the formation of such rings may, instead of being caused by 
strong bars, originate due to a weak oval distortion or the 
tidal effects of a companion (e.g. Combes 1988b, Buta et al. 
1995). Knapen et al. (2004) find a nuclear ring in NGC 278 
triggered by the non-axisymmetric potential induced by a 
recent minor merger. 

Finding dynamically induced starbursts in interacting 
galaxies is not unusual (e.g. Jogee et al. 2002; Keto, Ho 
& Lo 2005; ) where rings and inner spirals/bars are asso- 
ciated with the locations of resonances. A prominent and 
well studied example is the nearby starburst galaxy M 82 
which exhibits a molecular ring with a radius of 200 pc 
as seen in CO 1-0 (Walter, Weiss k Scoville 2002). This 
ring has also been studied in the highly excited CO 6-5 gas 
(Seaquist, Lee & Moriarty-Schieven 2006) where it appears 
to be dynamically associated with the transfer of gas from 
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x\ orbits in the large scale bar to perpendicular inner xi 
orbits. 



4.4.1. Evidence for a barred potential 

Work by Chapelon ct al. (1999) (C99) adresses the is- 
sues of starbursts in barred spiral galaxies. C99 claim that 
NGC 1614 is a barred galaxy with a bar PA of 37°, a bar 
length a of 5" and a deprojected a/b ratio of 0.23. C99 
based their conclusions on optical CCD images and it is 
interesting to note that this PA is similar to that for the 
crossing dust lane and CO distribution. 

Limits on the parameters of a possible bar in NGC 1614 
can be put on the onset of the prominent two-armed spi- 
ral structure. Defining where the arms actually start is not 
straight forward - but an estimate from the HST image sug- 
gests 7" as a projected limit. The prominent dust lane cross- 
ing just north of the central regions of NGC 1614 makes it 
difficult to judge the stellar structure and backbone of the 
galaxy in the visible band. We therefore investigated a 2.2 
micron K-band image with 0.3 arcsecond resolution (from 
Calzetti 1997). The K-band image reveals that between a 
radius of 1.1 and 0.4 kpc there is an apparent oval distortion 
of a/b axis ratio of 0.66 (and PA of degrees) dropping to 
unity at radius 0.4 kpc. This is roughly where the nuclear 
ring commences. There is thus either a weak north-south 
bar with projected length 2.2 kpc, or there is a warp in the 
inner kpc. Warps can occur in interacting systems. One ex- 
ample is the warped gas of NGC 3718 (Krips et al., 2005). 
One difference here is then that the dust lane of NGC 1614 
would be at an angle to the warp - and not part of the 
warped structure - as it is in NGC 3718. If there is a warp 
in NGC 1614 it would have to be gentle (because of the 
high axis ratio) but still happen rather quickly (from 1.1 
to 0.4 kpc). Thus, instead of a bar with PA 37° suggested 
by C99 we propose a north-south bar based on the K-band 
image. 

The CO velocity field suggests an S-shaped distortion in 
the iso-velocity contours. This is typical of barred galaxies 
- but may also be the results of a warp (or ouflows or other 
non-circular motions). The CO morphology is strongly af- 
fected by the crossing dust lane but is consistent with a 
north-south structure. It is important to note that the CO 
peak and distribution is offset from that of the starburst 
ring. The CO peak is located 1" north-west of the radio con- 
tinuum peak (and center of the radio continuum ring) . This 
is to some degree caused by the crossing dust structure - 
but not entirely. There is also a compact feature at w=4870 
kms _1 that helps to shift the integrated CO peak to the 
north-west. This could be the result of orbit crowding, but 
higher resolution observations are necessary to determine 
the inner CO dynamics of NGC 1614. It is furthermore in- 
teresting to note that in both the CO and Hi velocity field 
the velocity contours make a twist at the location of the 
starburst ring resulting in a local PA of 30-40 degrees (see 
Fig. [5]) . Such a PA is consistent with the shape of the Pa 
ring morphology. Again, high resolution CO data are neces- 
sary to determine both the CO morphology and dynamics 
across the ring and in the surrounding structure. 




Fig. 10. Overlay of the main dynamical features (repre- 
sented with arcs) over a Paa map (same as in Fig 5) in 
greyscale. The solid black arcs represent the main spiral 
arms. The dashed arcs represent the leading edges of the 
main bar. The white arcs represent the putative inner de- 
coupled bar. 



4.5. Nuclear radio continuum morphology and star formation 

Figure H] (lower right panel) shows the 1.4 GHz MERLIN 
in greyscale and the 5 GHz MERLIN in contours. The star 
forming ring is clearly seen in the 5 GHz data, while the 
nuclear region is better shown in the 1.4 GHz data, which 
has one of its peaks close to the center of the star forming 
ring. 

The strongest 1.4 GHz peak (the northern one) is 5.5 
mjy in a circular beam of 0."3. This can be converted to a 
source temperature of 53 000 K. This is an upper limit for 
the central peak, at this resolution, and is consistent with a 
scenario of a cluster of SNRs in the ring. However, an AGN 
can not be ruled out based on source temperature at this 
resolution, but higher resolution is needed. 

NGC 1614 has been observed with Very Long Baseline 
Interferometry (VLBI) by Corbett et al. (2002), who did 
not find any unresolved radio core. Hill et al. (2001) used 
long-baseline interferometry with The Australia Telescope 
and reported that less than 0.9% of the flux density at 2.3 
GHz is associated with an unresolved radio core. 

In Fig. [4] the 1.4 GHz continuum is more centrally 
peaked as compared to the 5 GHz continuum, which instead 
shows a ring with a radius of ~ 1". We interpret this as that 
the 1.4 GHz continuum is mostly related to supernova rem- 
nants, and traces old star formation, while the 5 GHz and 
8.4 GHz continuum is more related to free-free emission 
from HH-regions, and thereby traces current star forma- 
tion. The ages of the starbursts can be estimated using the 
spectral index, a, (assuming Soc v a ) as calculated between 
the 1.4 GHz and 5 GHz maps (the upper panels in Fig. [4]) 
toward the relevant areas (e.g. Hirashita & Hunt, 2006). 
Toward the center the spectral index is a ~ —0.9, which 
indicates clearly synchrotron dominated emission and an 
age of the starburst exceeding ~ 10 Myr. Although the fi- 
delity of the MERLIN 1.4 GHz map is poor, the spectral 
index in the star forming ring is significantly different and 
is a ~ 0, which indicates thermal emission and a starburst 
age of ^ 5 Myr. 
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We conclude that our radio continuum observations are 
in good agreement with the notion of an old nuclear star- 
burst, and younger starbursts in a R = 300 pc ring. 

4.6. Is there an AGN? 

We see no direct evidence of an AGN in our data. However, 
the presence of an AGN can not be ruled out, and is not 
in direct contradiction with our observations. The middle 
peak in the MERLIN 1.4 GHz (and in the MERLIN + 
VLA 5 GHz) radio continuum appears close to the center 
of the star forming ring, and could be associated with an 
AGN. Observations at higher angular resolution and sen- 
sitivity are needed to resolve this issue, although neither 
Hill et al.(2001) nor Corbett et al. (2002) found any clear 
indications of an AGN (see sect. I4.5[) . 

Risaliti et al. (2000) report on an X-ray nucleus 
obscured by a Compton thick column of gas (Njj > 
10 24 cm -2 ), based on the power law of the X-ray contin- 
uum. However, low signal to noise measurement of a power 
law continuum is difficult to interpret, and X-ray binaries 
have power law energy distributions that can mimic an 
AGN (e.g. Ebisawa et al., 2008). Wilson et al. (2008) have 
observed NGC 1614 with the submillimctcr array (SMA) 
and suggests a nuclear non thermal component based on ra- 
tio of submillimeter fluxes. However, the non thermal com- 
ponent could be associated with either an AGN or with 
SNRs. 

The central peak in our CO position-velocity diagram 
(Fig. right panel) could be associated with a nuclear, 
Compton thick region, obscuring an AGN. We have calcu- 
lated the largest possible radius of a Compton thick region 
(Nh > 10 24 cm~ 2 ) with a spherical geometry. We assumed 
that at most 30% of the flux (corresponding to 4.0 x 10 8 
Mq) in the position- velocity diagram (the central peak) 
could be associated with such a region. The largest possible 
radius, with N H > 10 24 cm" 2 would then be 140 pc (0."5), 
and it would not be resolved in in our CO data. Higher an- 
gular resolution would be needed to resolve the Compton 
thick region. However, based on the lack of obvious signs 
of an AGN, we conclude that the LINER-like optical spec- 
trum of NGC 1614 is probably related to nuclear starburst 
activity, and not to an AGN. 

Alonso-Herrero et al. (2000) has described a scenario 
where a LINER-like optical spectrum would be the con- 
sequence of the shock-heating by supernovae of the sur- 
rounding gas after a burst of star formation. This would 
typically occur 10-15 Myr after the onset of the starburst. 
NGC 1614 is likely an example of this scenario, where the 
nuclear (45 pc) starburst has been reported to have an age 
> 10 Myr (AH2001). The outflow, or super wind, described 
in Sect. 14.21 may also be related to the shock fronts giving 
rise to the LINER-like spectrum. A similar scenario was 
suggested to cause the LINER-like spectrum in the LINER 
galaxy NGC 5218 (Olsson et al., 2007). 

4.7. Feeding of the central region 

In the position velocity diagram of the CO-data (Fig. [3J 
right panel), three peaks are resolved. The southern and 
northern peaks (associated with the star forming ring) are 
stronger than the central peak. We estimate that 25% of the 
gas mass in the central beam (corresponding to 3.3 x 10 8 



Mg) is associated with the central peak. This suggests that 
there has been a substantial inflow of gas to smaller radii 
than the star forming ring at R = 1" in NGC 1614. 

Based on the large scale optical and Hi morphology, 
with a tidal tail and plumes, we speculate that the tidal 
interaction has triggered the formation of a bar where the 
gas has been transported along the leading edge and piled 
up in the star forming ring. An inner bar has decoupled dy- 
namically at this radius, and allowed the gas to flow toward 
the center. Star formation has occurred in this dynamically 
decoupled bar, and can be seen in the 1.4 GHz MERLIN 
map (Fig. [3J . Inflow occuring before the formation of the 
ring, and the continued inflow along a decoupled bar af- 
ter the ring formation, has triggered and maintained the 
nuclear activity of NGC 1614. 

4.8. The nature of the crossing dust lane and the interaction 

Prominent dust lanes in merging galaxies are not uncom- 
mon and have a variety of properties. One example is 
the dust lane of the NGC 4194 merger (e.g. Aalto and 
Hiittemeister 2000). Just like the NGC 1614 dust lane, it 
is a minor axis dust lane crossing the galaxy at an angle to 
the velocity field. 

In the NGC 4194 merger it is speculated that the dust 
lane is feeding the central region with H2 for star formation 
from the outer regions of the merger. A similar scenario is 
possible for NGC 1614 and should be studied at higher 
resolution. It is interesting to speculate whether the minor 
axis dust lane is part of an inner polar ring structure. 

In a numerical simulation study of minor- or interme- 
diate mergers, Bournaud et al. (2005) found that the gas 
brought in by the disturbing companion galaxy is generally 
found at large radii in the merger remnant. The gas is re- 
turning to the system from tidal tails and often forms rings 
- polar, inclined or equatorial - that will appear as dust 
lanes when seen edge-on. 

The CO emission we identify as associated with the fore- 
ground dust lane is blueshifted with respect to the expected 
rotational velocity at the same position in the plane of the 
galaxy. This is consistent with it not being in the same 
plane, but does not constitute evidence. Alternatively the 
dust lane acts as a gaseous bar, inclined to the K-band 
stellar bar, along which gas is being funneled to the center. 
This interpretation is not consistent with the CO velocity 
field however. Higher resolution molecular observations are 
necessary to distinguish between the polar ring/gas bar sce- 
nario and to investigate how the gas is flowing towards the 
nucleus. 

Both a polar ring and a bar can be triggered by the 
interaction just like the tidal features of NGC 1614. Neff et 
al. (1990) suggest that an intruding companion is located 
10" to the south-west of the center of NGC 1614. Peculiar 
Ha velocities and an alignment of the jet-like optical tail 
suggests this. An inspection of the HST WFPC2 F606W 
image does indeed reveal a drawn-out feature here pointing 
in the direction of the tail. We find no molecular counter- 
part to this object, but this is hardly surprising if it has 
already been stripped of its gas content and only has its 
nucleus left. The feature also appears in the K-band image. 
AH2001 suggest that there is evidence of the nucleus of a 
companion approximately 1" north-cast of the nucleus in 
their B-band image. There is no hint of this companion in 
the K-band image - and no evidence of a particular dis- 
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turbance in our CO or HI velocity field here, although the 
spatial resolution needs to be higher to be certain. 

5. Conclusions 

1. We detect 3 x 10 9 M Q of molecular gas within a north 
to south slightly elongated bar like structure with a size 
of - 2.2 x 1.5 kpc. 

2. The north to south rotation of the central kpc molec- 
ular gas is consistent with an inclined disc / bar that 
follows the expected rotational direction based on the 
assumption of trailing spiral arms. 

3. There is a molecular extension to the northeast, associ- 
ated with a crossing dust lane. 

4. The molecular gas in the central kpc is double peaked, 
with peaks at R = 300 pc. There is a third, central peak 
at lower integrated intensities, probably associated with 
a nuclear starburst/AGN. 

5. We have detected a MERLIN and VLA 5 GHz radio 
continuum ring with a radius of 300 pc. This ring is 
coincident with previous radio continuum and Paa ob- 
servations. We conclude that the radio continuum ring 
originates in recent star formation. The CO peak inten- 
sity is shifted 1" to the north-west of the center of the 
starburst ring. 

6. The MERLIN 1.4 GHz radio continuum is triple peaked, 
with a peak separation of ~ 1" (300 pc). The brightest 
peaks occur in the northern and southern part of the 
star forming ring with a radius of 300 pc. The other 
peak, which is almost as bright, occur inside the star 
forming ring, close to its center. The brightness temper- 
ature of maximum 53 000 K of the peaks are consistent 
with a cluster of supernova remnants, although an AGN 
can not be ruled out for the peak close to the center. 

7. The position velocity diagram of the Hi-absorption was 
used to calculate a dynamical mass of 3.5 x 10 9 M in the 
central R — 1" (300 pc). An upper limit to the molecular 
gas mass in the same region is ~ 1.3 x 10 9 M©. This gives 
a gas mass fraction of <~ 0.4 or less in the central R = 1" 
(300 pc). 

8. A solid body rotation curve was fitted to the Hi- 
absorption velocity field, out to a radius of 2" (600 pc). 
The residuals are blueshifted over all of the fitted area. 
We interpret this as an outflow or a superwind, with a 
deprojected velocity of ~160 kms _ . 

9. We suggest that gas may flow to the nucleus along a 
dynamically decoupled bar inside the star forming ring 
(or spiral), which we see tentative evidence for in the 
1.4 GHz MERLIN map. The existence of such a bar 
should be investigated with CO observations and higher 
angular resolution. 

10. The LINER like spectrum is likely related to the shock 
fronts associated with either the large scale super wind, 
or with shock fronts on smaller scales related to su- 
pernova remnants from previous bursts of star forma- 
tion. We conclude that the LINER activity observed in 
NGC 1614 is probably due to starburst activity, and not 
to AGN-activity. 
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